Shot noise in superconducting junctions with weak link formed by Anderson impurity 
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A theory is developed to study shot noise in superconducting (SAS) and hybrid (SAN) junctions 
with singly occupied Anderson impurity (A) as a weak link. The zero-frequency dc component of 
the shot noise spectral density is calculated at zero temperature as a function of the bias at different 
Coulomb repulsion strengths ([/), and show a remarkable structure resulting from combination of 
electron-electron interaction and Andreev reflections. 



Motivation and main results: Following the first the- 
oretical work on noise in superconducting point junctions 
the underlying physics attracted a considerable the- 
oretical Ipj- fzl an d experimental attention (see re- 
view article for more complete list of references) . In 
most of these works, the quantum point contacts have 
been considered in both the linear and nonlinear regimes. 
Tunneling junctions in which the electrodes are sepa- 
rated by a resonant double barrier were studied in Refs. 
|^,|^. Novel experimental techniques now enable the 
study transport and shot noise in quantum dots located 
between either normal (N) or superconducting (S) elec- 
trodes. So far, however, the physics of shot noise in su- 
perconducting junctions with strong electron correlations 
has not been exposed. Investigation of this fundamental 
aspect is carried out in this work. The shot noise spectral 
density in SAS and SAN junctions is calculated below 
for the case were the Anderson type impurity level {A) 
is singly occupied. For an SAN junction at small bias 
voltages V, doubling of the normal Poisson noise to cur- 
rent ratio (Fano-factor) is preserved, although its depen- 
dence on the electron-electron interaction is quite essen- 
tial. In an SAS junction, the main process contributing 
to the current and shot noise power is multiple Andreev 
reflections (MAR). However, in spite of the fact that 
for low V the number n w 2A/eV of MAR is large (A 
is the superconductor gap), the current and shot noise 
density are rather weak. This is due to the low effec- 
tive transparency F of the junction as a consequence of 
interaction (Coulomb blockade). Large-n processes are 
therefore damped as F". 

Model and eflfective action: For convenience, the 
junction is represented by two half electrodes on the left 
(L) and the right (i?) separated by quantum dot. The dot 
(located at the origin) is modeled as an Anderson impu- 
rity A with level position eq < and Hubbard repulsion U 
under the condition for single occupancy t/ > — eq > 0. 
This can be justified by noticing that in recent exper- 
iments pj|j20[| on semiconductor quantum dots it was 
shown that tunneling takes place through a separate state 
with the sign of a Kondo behavior (the tunable Kondo 
effect). 

The starting point is the tunnel Hamiltonian H — 



Hl + Hn + Hd + Ht in which Hj{i = L,R) are lead 
Hamiltonians (usually of BCS form) defined in terms of 
electron field operators tpj„{x,y,t), = co Scr '^cr'^o- + 
Un^ni and Ht = ^X^jo- cJ.V'jCT(0, i) -I- he. Formally, 
the whole physics is contained in the partition function 
Z = J d[F]exp(iS) where the path integral Grassman in- 
tegration is carried out over all fermion fields [F] and the 
action S is obtained by integrating the Lagrangian per- 
taining to the Hamiltonian H along a Keldysh contour. 
In a recent work |l5) we studied I — V characteristics in 
SAS and SAN junctions, and developed a formalism to 
carry out integration over the fields ipja {x, y, t) , combined 
with the dynamical mean field approximation to account 
for the quartic term Un^ni. At the end of this procedure 
one arrives at an effective action Sef / which depends on 
a parameter 7_ whose physical meaning is the difference 
between spin up and spin down energies of the dot, (to 
be determined self consistently in terms of the effective 
action). The latter is defined in terms of lead Green 
functions, gf^^^^ where j — L. R and R/A/K denotes 
advanced, retarded and Keldysh respectively. Explicitly, 
dropping the lead index. 



is a matrix in Keldysh space with 



{e±iO)T, + i\A\Ty 



(1) 



(2) 



g^(e) = (.9^(6)- 5^(6)) tanh(e/2T). (3) 

The set of Pauli matrices Tx,y^z acts in spin (Nambu) 
space, and the lead dependence enters through the cor- 
responding superconducting gaps A. It is now possible 
to define the kernel, 

G"i(e, e') - S{e - e'){e + 7_ - r,?) + zTg+{e, e')T,, (4) 

(where g+ = {g^ + 5fl)/2) as an operator in Nambu 
Keldysh product spaces (here e = eg -t- [//2 and F oc is 
the usual transparency parameter). The effective action 
is a functional of the Grassman fields c, c which are four 
tuples in this space. 
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S'off = y — J de'cG ^(e,e')c. 
The self-consistency equation for 7_ is, 
7_ — < / dtca^c >, 



(5) 



(6) 



where <^ O 



and the set of Pauh 



J" d[c]Oexp{iS f) 
J d[c]exp{iSef f) 

matrices <Jx,y,z acts in Keldysh space. 

The physical justification for using the mean field ap- 
proximation is related to the fact that the the Kondo 
resonance is suppressed by the superconducting order pa- 
rameter. Such suppression occurs due to strong attenua- 
tion of the density of states in an energy region of order 
A near the Fermi energy. Thus the number of low energy 
electrons which are able to screen the local impurity spin 
is small p6| , pT[ . Therefore, we can expect that ineffective 
screening results in a larger domain in parameter space 
(eo, U) for which the single occupancy (doublet state) be- 
comes the ground state of the whole system. This state 
is obtained from the solution of the mean field equation 
(i). 

Explicit expression for shot noise; The noise spec- 
trum measures the current fluctuations in the junction. 
It is defined by the symmetrized current-current corre- 
lation function which, in terms of current operators in 
Nambu (8) Keldysh space reads, 



+ <f/(i)(i2)/('>(ii) > -2</>2], 



(7) 



where T is the time ordering operator and < ... > denotes 
quantum mechanical thermodynamic averaging with re- 
spect to the total Hamiltonian H. Starting from the 
general definition of the current operator for tunneling 
through a quantum dot [ISl, in the present case it reads. 



/(1,2) 



±- 



En 



ijRk{0) - h.c], 



(8) 



where ipRki^^) is the Fourier transform of 7/'fl(0,y,t) (the 
so called surface field operator) with respect to y. Substi- 
tuting the explicit form for the current operators ^ into 
equation (|^) we then obtain an expression for K{ti,t2) 
which involves Grassman integration over surface fields 
and dot electron operators. The first integration is Gaus- 
sian and can be done exactly, involving the Green func- 
tion matrix g. Integrations over the dot fcrmion fields is 
completed within the dynamic mean field approximation. 
After Fourier transform on ti — 12 it is possible to express 
the power spectrum A'(w) and 7_ in terms of the Green 
functions of the entire system G^'^ (inverse of the ker- 
nel defined in Eq.(^). In the present work, attention is 
focused on the zero-frequency dc component of the shot 
noise K = K(0) at zero temperature. The first novel 



result of our study then consists of a workable expres- 
sion for the noise power density function K = Ki -\- K2 
supported by a self-consistent equation for the energy 
occupancy parameter 7_, that is. 
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g^){G'^^G-)~g^F]}, (9) 
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where, 



F 



-iV 
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(11) 



(12) 



(13) 



(14) 



We arrive now at our main goal, namely, using expres- 
sions (p[]l^) to analyze the noise in SAN and SAS 
junctions. Before starting this analysis we note that 
Eqs. (§,^ 11) can be easily generahzed for the case 
with different left and right transparency parameters: 
Tji,Tl. To do this we write Eqs. (|s|JlO|) in a sym- 
metric form. Moreover, Eq. (O) is to be replaced by 
F = ^G^iTng^ + TLg^)T,G^ 

Shot noise in SAN junctions: In this case only one 
Andreev reflection takes place and K2 can be written 
explicitly as an integral over the whole energy domain 
(including sub-gap), while for JsTi, integration is effected 
above A only. From that part of Ki + K2 related to in- 
tegration with |e| > A we easily obtain the shot noise for 
NAN junction, (tunneling between normal metal elec- 
trodes through an Anderson-impurity center). It assumes 
a standard form [Q, although it depends on the interac- 
tion through the parameter 7„, that is. 



K = 



[r_(i-T_) + r+(i-T+)], 



(15) 



where T± = r-^/((e ± 7_)'^ -I- F'^). In a perfect resonance 
condition (e = and ?7 = 0) we get K — as expected 
for a pure point junction iQ,^. When U ^ we have 
to solve the self-consistent equation (11) to find 7_ and 
substitute the solution for a given voltage into equations 

Another limiting situation for which a simple analyt- 
ical expression exists is that of small bias eV < A. For 
completely transparent point junctions the current-noise 
spectral density vanishes m . For non-resonant transport 
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there is a finite shot noise in the sub-gap region which, 
in the hmit of low effective transparency, yields the Fano 
factor to be equal to 2. The principal contribution to 
the noise density in the regime eV <C A comes from K2. 
Therefore, integrating over the sub-gap energies in equa- 
tions |ll|) results in, 

2e^\V\ r4[r2?2 + (^2 _~2)2] 



K 



TTh [r27i + (g2 +0.5r2 -72)2]2- 



(16) 



This formula, together with the self-consistency equation 
( |rT| ) yields the shot noise density for T = and eV <^ A. 
For r ^ e, 7_ it leads to twice the normal Poisson noise- 
current ratio. Numerical calculation using the exact ex- 
pressions (|(lO) for specific values of U indicates that 



the deviation from the approximate relation (16)) starts 
already at eV ^ A/2. Therefore, even in the sub-gap 
regime, the non-linear ^-dependence of the shot noise is 
important and exact expressions for the spectral density 
should be employed. The shot noise of an SAN junc- 
tion is displayed in Fig.l (hereafter, all the energies are 
expressed in units of A ). 
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(17) 



where N(e) - ^j,^. 

Shot noise in SAS junctions: The current-noise spec- 
tral density of an SAS junction is calculated numerically 
from equations ([9HlO|). For a constant bias V it is use- 
ful to discretize the energy integration JTsf . The energy 
domain in equations (|9|, p^ , [rT| ) is divided into slices of 
width 2eV and integration is performed on an interval 
[0 < E < 2eV]. The Green functions become matrices 
with indices representing the different energy slices. Ex- 
plicitly, 



{m\G-\e)\n) ^ S„ 



fm + 7- - Tze+ —gnX'^m )Tz 



y("^l5L(e)|n)r^, 



(18) 
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FIG. 1. The shot noise K versus the bias V for an SAN 
junction at sub-gap voltages. The parameters are eo = —1.5, 
r = 0.6, U = 2.45 and 2.713. 

For the parameters employed in the figure, the noise in 
the sub-gap regime is rather small, starting to grow near 
the superconducting gap. It is strongly dependent on 
the Coulomb interaction which, as a general rule, sup- 
presses it in comparison with its value appropriate for 
simple tunneling through a non-interacting impurity. It 
is interesting to note that the role of the repulsive Hub- 
bard constant U in an SAN junction (in the regime of 
single occupancy) is similar to that of an exchange term 
in a contact between a superconductor and a ferromag- 
netic metal. Therefore, our results might be relevant for 
this junction as well. Note also that for clean point junc- 
tions, there is a saturation value for the shot noise spec- 
trum density at large voltages [Q. This value is read- 
ily reproduced from our formulas (|^,^) if we consider 
the limit F A at resonance (e = 0) for bias volt- 
age eV > A. The only contribution to the correlation 
function comes from energy integration above the gap in 
equations ([9|pO|), yielding the saturated value, 



9l (em)T+ (19) 

where e^n — e + 2meV, = ± eV, the super- 
scripts denote the matrix elements in Nambu space and 
P± = (1 ± r2)/2. We dropped the constant phase differ- 
ence (f)o = arg(Ai — A/j) which does not play a special 
role here. 

We are mainly interested in the voltage dependence at 
the sub-gap region where MAR are important and the 
way how this noise- voltage characteristic is influenced by 
Coulomb interaction. The result is displayed in Fig. 2. 




FIG. 2. The shot noise K versus bias V for an SAS junc- 
tion at sub-gap voltages. The parameters are eo ~ —1.5, 
r = 0.6, U = 2.4 and 2.7. 

As can be deduced, there is a strong effect of the 
Coulomb repulsion on the noise spectrum K in the sub- 
gap region. For voltages not too close to 2A the noise is 
smaller for the higher value of U. This is similar to what 
happens with the current itself, as can be judged by a 
glance at the I — V curve (Fig. 3). 
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FIG. 3. I — V curves for an SAS junction at sub-gap volt- 
ages. Parameters are the same as in Fig. 2 

And yet, the Fano factor K/2eI is nearly independent 
on interaction (see Fig. 4). This finding stresses the im- 
portance of MAR process, analogous to what has been 
argued in junctions with low transparency |q]. 




FIG. 4. The Fano factor K/2eI versus bias V for an SAS 
junction at sub-gap voltage. The parameters are as in Fig. 2 

As in the case of SAN junction there is a saturated 
value for the spectral noise density power also for SAS 
junctions. In the limit F 3> A this saturated value is 
twice of that of Eq.(0). 

In conclusion, based on a theory developed for the 
study of tunnehng in SAS and SAN junctions, the 
shot noise in such systems is calculated. Special attention 
is devoted to analyzing the implications of the Coulomb 
repulsion between electrons in the dot on the tunneling 
process in general and the noise spectrum in particular. 
The theoretical treatment uses a combination of Keldysh 
non-equilibrium Green function and path integral formal- 
ism, and the interaction is handled within the dynami- 
cal mean field approximation. General expressions for 
the current noise spectral density correlation function are 
then derived. 

The main results of the present research can be sum- 
marized as follows: 1). In SAN and SAS junctions, the 
Coulomb interaction results in a nonzero value for the 
occupancy energy (7_) and thus always acts as a factor 
which lowers the transparency. Therefore in SAN junc- 



tions the shot noise appears already at eV < A (even if 
e = 0). The doubled value of the Fano factor emerges in 
the limit F ^ e, 7_. 2) In SAS junctions, the shot noise 
correlation function in the sub-gap regime {eV < 2A) is 
suppressed at higher values of U. The Fano coefficient is 
nearly interaction independent, and unambiguously man- 
ifests the importance of M AR. 
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